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Use of a cDNA microarray to 
analyse gene expression 
patterns in human cancer 

Joseph DeRisi**, Loliia Pcnland- & 
Patrick O. Brown- (Group 1); 
Michael L Bittner*', Paul S. Melt2er\ 
Michael Rav-^. Yidong Chcn\ Van A. Su-^ & 
Jeffrey M. frcnf^ (Group 2) 

The development and progression of cancer^"^ and 
the experimental reversal of tumorigenicit/-^ are 
accompanied by complex changes in patterns of 
gene expression. Microarrays of cDNA provide a 
powferful tool for studying these complex phenom- 
ena^®. The tumorigenic properties of a human 
melanoma cell line, UACC-903. can be suppressed 
by introduction of a normal human chromosome 6, 
resulting In a reduction of growth rate, restoration 
of contact inhibition, and suppression of both soft 
agar clonogenicity and tumorigenicity in nude 
mice^-^ ®. We used a high density microarray of 
1,161 DNA elements to search for differences in 
gene expression associated with tumour suppres- 
sion in this system. Fluorescent probes for 
hybridization were derived from two sources of cel- 
lular mRNA [UACC-903 and UACC-903(+6)] which 
were labelled with different fluors to provide a direct 
and internally controlled comparison of the mRNA 
levels corresponding to each arrayed gene. The flu- 
orescence signals representing hybridization to 
each arrayed gene were analysed to determine the 
relative abundance in the two samples of mRNAs 
corresponding to each gene. Previously unrecog- 
nized alterations in the expression of specific genes 
provide leads for further investigation of the genet- 
ic basis of the tumorigenic phenotype of these cells. 

DNA micnnirrays, conuunin^ 1,161 tola! elements, 
includini; S70 different cDNAsiUui ccnilrols**"" (sec 
Methods), were printed roboiically onto a i;lass micro- 
scope slide in lour quadrants coverini; an area of about 
I cm- (Fig. I ). We prepared tluore.scent cDNA probe.s 
using total pol)' ( A)* mRNA from UACC-9()3 cells and 
UACC-903(+6) cells by labelling with a green and red 
fluor respectively. A mixture of the two flourescently 
labelled probes was hybridi/.ed to the DNA inicroarray. 
This comparative hybridisation method, coupled with 
the doping of synthetic standards and an cstiniation of 
statistically significant deviation for local background 
variance allowed a direct and quantitative comparison 
of the relative abundance of individual PNA sequences 
in this complex sample'*"**. We added a set of synthetic 
poly (A)*-tailed 'niRNAs* to the purified mUNA from 
each cell line as internal standards to assist in quantita- 
tion and estimation of experimental variation intro- 
duced during labelling and reading. Targets 
complementary to these standards were included, in 
duplicate, on the micmarray. Based i>n these standards. 
mUNA species comprising 1:10,000 of the mass of the 
poly (A)* UNA could readily be detected. 

In a representative two-colour fluorescent scan of all 



spond to genes preferentially expressed in the tumori- 
genic UACC-903 cell line, and the reddish spots corre- 
spond to genes preferentially expressed in the 
non-iumorigenic UACC-903(+6) cell line. Genes 
expressed at approximately equal levels in the two cell 
lines appear yellow or brown. A portion of the array at 
higher magnification highlights the diverse pattern of 
differential expression observed (Fig. 2b). In Fig. 2c, rec- 
tangles corresponding to specific array elements are 
coloured to reproduce the hue and intensity of the fluo- 
rescent signal at each element. The hybridization signals 
from a duplicated set of genes are shown juxtaposed, to 
illustrate the reproducibility of the hybridization signals 
for each gene. 

To address the possibility that an apparent difference 
in expression might result from experimental variables 
unrelated to the difference in chromosomal composi- 
tion between the two cell lines, we examined the vari- 
ance in expression for 90 ^housekeeping' genes. We 
selected these genes based on the assumption that they 
would not be differentially expressed between the two 
cell lines. The averaged red/green ratio for this subset of 
genes was 1.13. The averaged red/green ratio for the set 
of five internal standards was 0.97 (ii = 10). The vari- 
ability in the expression level of the housekeeping genes 
probably overestimates the experimental variability in 
measuring differential expression. As a conservative stan* 
dard, an absolute fluorescent signal (red or green) with 
an intensity greater than that observed at the control 
array elements containing total human genomic DNA 
was considered to represent specific hybridiz;ition. Gene- 
specific hybridization was therefore only considered sig- 
nificantly different beiween ,s;unplcs if the following two 
criteria were met: i) the signal intensity (green or red) 
exceeded this threshold; and ii) the logarithm of the 
red/green lluorescence signal ratio differed by >3 S.D. 
from the nK*an logarithm of this ratio for the *house- 
keeping* gene panel (that is, ratios <0.52 or >2.4). 

By these criteria. mRNA levels lor 1 5/R70 (1.7%) genes 
were significantly diminished, while the mRNA levels 
for 63/S7() (7.3%) i:enes were significantly increased in 
association with suppression of tumorigenicity by intro- 
duction of chromosome 6. To test the reliability of 
microarray hybrid i/.;U ion results in ideruifying differen- 
tially expressed genes, we analysed 16 genes by north- 
ern analysis. In each case, the results of northern analysis 
corroborated the differential gene expression identified 
by microarray hybridization (Fig. 3). 

Significant differences in expression between these 
two cell lines identified several genes as candidates for 
determining features of the tumorigenic phenotype of 
the melanoma cells. For example, among the genes 
detected with significantly higher expression (> 10-fold) 
in the tumorigenic cells was the human brown locus pro- 
tein (TRPI/melanoma antigen gp75). This is the most 
abundant glycoprotein in mclanocytic cells and a critical 
melano.some membrane protein '-*' \ Additionally, its 
expression is reduced when melanoma cell lines arc 
induced to differentiate by treatment with HMBA^^''\ 
Also expressed at a significantly higher level was a spliced 
variant of the mRNA encoding myelin P!.P/l)M20.This 
is widely expressed in neural crest derived cells in early 
development and has been suggested to play a role in 
cell-cell signaling during development'"*. 
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Rq. 1 Properties of cDNA nwcroarrays. «. A fluorescent scan of DMA printed onto a po»y.|ystne coated sUde. The DMA is stained with a DNA-soecific ftuorescent 
dye. YOYO. The cemer-to-center spacmg of adjacent spots is 450 m. altowtng the potential for up to 10.000 spois/2,54 X 7.62 cm microscope slide b Effi- 
cient blocking of hybridization to DNA repeats. Hybridization of fluorescein- labelled poly (dT)* to arrays in the absence of competitor produces strono 
hybridization to immobilized poly (dA)* as well as to some cDNAs. such as the EST T64827 shown. Rhooamine-iabeited cDNA (red) from the UACC-903 c^ 
tine hybndized m the presence of poly (dAp blocker shows irttle if any signal at either site fTotai H = total human). StmUarty. hybndization with fiuorescem-labelled 
Coti DNA in the absence of competitor produces bright signal on immobiltzed Coti DNA, total human DNA and at some cDNA elements (presun^ed to con- 
tain highly repeated sequences, such as R23416): while Rhodamtne-iabelled cDNA (red) from the UACC-903 cell line produces little it any stonai at these 
locations when hybndized in the presence of excess unlabelled poly (dA)- . and human Coti DNA The absence of signal at some cDNA locations foliowtnc 
UACC-903 cONA hybndizations also indicates that the PCR-amplified. ptasmid vector sequences at all cDNA tatgets oo not contribute siqmlcant hybndiz^ 
tion signal, c. Schematic of the array organisation. Robotic printing from 96 well microtiier trays was earned out with 4 prmt heads spaced to fit into 4 adia- 
cent microtiter wells. This maps the contents of each tray into four separate quadrants on the glass slide. A colour-coded map of the general distribution of 
target types in each of the resulting quadrants IS shown. 



cIs were clcviiicil by the ;iddilit>n i>l a norniitl chrnmo- 
somc 6(17 i;cncs) arc kmnvn lo be activated by IFN-y. a 
cardinal proinllainmatory cytokine that, anitinj; other 
activities* induces e.\prt*ssion of the gene products of the 
MHC class II locu.s. For example, the mUNA encodinj; 
monocyte chemotaclic protein I (MCAF/MCPl),a 
cytokine that induces monocyte chemoiaxis and activa- 
tion''* '^ was more than 10-fold less abundant in the 
lumorigenic cell line. In the skin. MPCl is critical in the 
regulation ofcuianeous monocyte traffickini; and 
elevated expression plays a role in suppression of tumour 
growth and metastasis'''"-'. The mechanism by which 
these intert*eron-7 regulated genes are induced in UACC- 
903 cells by transfer of a normal chromosome 6 remains 
to be determined. It is worth noting, however, that the 
interferon-7 receptor gene is localized to the distal long 
arm of human chromosome 6. 

Finally, several genes that showed > 1 0-fold higher 
expression in the suppressed UACC-903(+6) cells have 
previously been recognized in other models of tumour 
suppression. Most notably, there was elevated expres- 
sion of the mRNA encoding WAFl (p2l ), a key media- 
tor of tumour suppression by p53 (rcf. 18). The p21 
protein had previously been identified as a melanoma 
differentiation-associated antigen (termed mda-b)''' -^. 
In melanoma cell lines suppressed for metastasis by the 
introduction of chromosome 6, e.xpression of WAFl 
(p21 ) mRNA and protein correlates inversely with 
metastatic potential-". 



These results provide a wide view of the diverse sys- 
tems that are altered in this model system of lumori- 
genicity, and focus attention on specific gene products 
and pathways that may he of particular importance in 
this tumour type. 

Our ability to classify human cancers in a way that 
reflects the underlying molecular pathology or that 
anticipates their potential for progression or response 
to treaimeni, remnins primitive. Using cDNA microar- 
rays to define alterations in gene expression associated 
with a specific cancer may be an efficient way to uncov- 
er clues lo the specific molecular derangements that con- 
tribute to its pathogenesis and thus identify potential 
targets for therapeutic intervention. Moreover, recogni- 
tion of pathognomonic alterations in gene expression 
might provide a basis for improved diagnosis and mol- 
ecular classification of cancers and thus allow selection of 
the most appropriate therapeutic strategies. 

Public databases of human expressed gene sequences 
contain partial sequences of at least 40,000 different 
human genes", and efforts to develop a human tran- 
script map have developed rapidly-'. Based on the high 
yield of information obtained using an array of < 1,000 
different genes, a more comprehensive survey of gene 
expression patterns, using a more complete array of 
human genes, will likely provide a rich source of new 
and useful insights into human biology and a deeper 
understanding of the gene pathways involved in the 
pathogenesis of cancer and other di.seases. 



^ CD 



letters 



o o o a a 

o n 0 

o o o .1 o .1 o 
i> *} o o J , 

o o o o - o 

o o o Q e 

o a J o 

3 



0 0 0 0 

0 0 
O 0 0 0 
0 0 0 0 

0 o o 
O O O 0 o 

0 • o 

0 


0 0 0 0 

0 0 o 
ODD 

-> 0 
0 0 

0 0 






5 


El 


- 


- 0 J 


1 nil 


IIIIIQI 


IIIBIII 



Fig. 2 DNA mtcroarray anatysis of changes in gene expression between the tumorigenic cell line. UACC-903. and its non-tumohgenic denvative. UACC903(4^. 
denved by introduction of a normal chromosome 6. a. A ratio image of the results of simultaneous hybridization of RhodaminellO-labelled cDNA (green) 
from UACC-903 and Cy3-labei»ed cDNA (orange-red) from UACC-903(-*-6) to a microarray. To oroduce this image, the scan images con-esponding to each flu* 
orescent probe were combined as the appropriate colour channels in a smgte image. Arrows indicate the location within the array of the corresponding 
genes analysed by nortr>ern btotttng (Fig. 3). b. A magnified image of the area of the array boxed tn white in (a), c. Magnified image of three cONAs indenti- 
fied by arrows tn (a), representing the cDNAs for: left. MCAF/MCP-I (r/g ratio >10): centre, f^acttn (r/g ratio 1.04); and right, n-l-antichymotrypsin (r/g ratio 
0.2) (see Fig. 3j. d. simptifted representation of ratio hybridization results. Quantitative fluorescence intensity data is extracted from each array target. Ttm aver* 
age target colour ratio determtr^es the hue of each box and the average intensity determines the bnghtness of each box. in this image, the order of the boxes 
corresponds to thetr original order in the microliter plate from which they were printed. Duplicate printings of the same plate can be examined side by side, 
as in the first two rows shown here, to assess reproducibility of the hybndization results (see text). Numbered arrows indicate the location withtn the array cor- 
responding to genes analysed by northern blotting in Fig. 3. 



Methods 

Generation of micruarrnys, hybridisation, scanning. 'Hk* 
prcp.u.iluin o( cn.ual inicri».s».ii|K .sliJcN .iiul Mili>a]UL*iu n»in»i- 
k priniini;nr |)NA w.is cjrricii tut! in ;i ni.iniui Miiiiki! i«i ih.ii 
ilocrilvit . HricOy, prc-ckMiial Nliik*> wwc irc.ttal with 
piiK -I -IvMiK' M»luiiiin (Stiini.ii tit lorm an .ulhoiw mm I.kv li»r 
priniiiij;. VCM prtichKlN. piirifial hy cllunol pnriric.iiioii. unv 
rcMispcniloii in 3.\ SSC A L'u>toni built iirrayint; rohm pivkal 
up anil ilcpo.Micil .small vtilumc> iijnoliior>) ol DNA onin 
ihc siulc>. Alter pnntinj:. the sliiic.s vwrc waviu'J in a I).2**/m .SDS 
.MiUiiMMi. *i he rcnuining hoiiiiif DNA \v.)> di-naiiitci) h> miIv- 
mcrjiini: the sliile.v in 95 ihstilieil water (*ir 2 tniti Inllnweil t>y 
a iTriel vvash with 9Vi„ elhanol. DNA was iioNslinked in the 
sliiic.s (Siratai:eno Siratalinker. ftO ml), lo prevent nnii-vt>ecirie 
pn»he hiiiiiinj:, the siiile> were hU»ekeif In rinvini: in a M»luiMin 
01*70 niM MiCLinie anhviiriiie ilissolveJ in D.I M hnrk aciii pH 
S.fl. ciintainini; I -nu'thyl-2-pvrritliiiinone (Atifnch). 

AtMitiiinal pnittKoU and parl> li>l pertaininj; l»» inieroarray 
taliricatitin can ho olitaincd fnini http://emj;fU.stanforii.edu/ 
phriHvn. 

Puriticd. ialvlled eDNA wa.s resuNpended in 1 1 p) ol .VSx SSC 
i(iniaiinni:4 (.ijiiil p«»ly (dA)* DNA. ^..'i h, t<i/MUNA,4 Mj;*)! 
huiiun CtM I DNA ((iiivd HUl.). and 1)..^ pi nl Mr.. SDV Prior to 
hyhridi/atiiin. the Miluiion ua\ hniled loi 2 niin then alhmwl to 
o»ol lo loom lemperalure. Ilyhi tdi/alion was carried out al 
b2 *C* lor * U h in a water hatli. I'rior io seannini:, .\lide> vveie 
washed in > SM f).y'„ SI >S lor 5 niin .uitl «).:>■ S^^C lor 1 inin. 

lot.al laser ninrtts^ope Iniilt l>\ S. .^nuth with >ottware wniien 



In N. /i\. A >epai aie N*;an. umui: the appropriate cuitatioi; line, 
was done loi eaeli ot the tno lUittrophores used. Data was col- 
lected .11 a inaxinuMi) lesoltition i>l S) niK'ittns/pt.vel with 12 bits 
o| ilepiii 

Probe preparation and laheliinj;. UNA was exiracied iVoni cclU 
UMiiv: ihelria/ol reauent (I'l l Inc. I. foUowinv: the manufactur- 
er > diivetions. i DNA probes were synihcsi/ed trom singly olipo 
li l - selected ( Piiariuacia ) inKN.A pooU. KUtttreMJcntly laMled 
el )NA was prepared Inuii inKNA by olij;o d l -priined |Milymcr- 
i/ation usinv: .SuperSertpt 11 reverse transcriptase (LTI Inc.). 
l he p«*ol ot nucleotuie.s in the lai>ellinj; reaction was 0.5 mM 
dC.TI'. dA'll' and dCVIP and 0.: niM dTTI». Fluorescent 
nucleotides. Khodaniine I H> dUTP (IVrkin F.limT Cetus) or 
C*.\ 3 dl'TP ( Aiuershani), were present at 0.1 mM. IVoIks were 
purified In i:el chroniaiojjraphy (HioSpin 6/litoKjd) and 
eihanol precipitation. 

Select inn of cDNA elements and generation of control tem- 
plates. Synthetic cl^NAs were prepared by clitning random 
IUiittH \ and Hi/idlll ended lrai:ntenis of F.. toli DNA in the vec- 
tor pSPM poiy (A)* (Prome\:al. lineari/iiig isolated plasmid 
DNA with K(i»KI and synlhesi/in^ p»»ly (A)* tailed UNA com- 
plenuniary to the insert lioin tlw resident SPO promoter 
(Pioineua). Prior to use. the s\nihesi/ed KNAs %vcre sclccied on 
olit*o d I' celIulo.se. The largest uioup of cDNAs conM.slcd of 674 
cDVA clones lioiii the I NIK anawd nonnali/ed infant brain 

i...f.i. . . I i»w'>c *.t.i:.. - ...... >;..wvt».J ti» i; » ,veiy Nli> 

hbrar\ nteinber that v«m tesj^oiuled to a named gene according 
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Ftg. 3 Nonnern hybridizaiion ftut>stamiai* 
tng the consisiency of the cONA mtcroar- 
ray results. Corresponding locations withm 
the cONA microarray illustrated tn Fig. 2a 
are provided lor VWat-l/p2l: 2) MARCKS: 
3) coWagenase: 4) MCAF/MCP-v, 5) n-l- 
anttchymotrypsin: and 6) -aatn. The sig- 
nal detected by a radio-iapeiled ^actm 
probe represents a control for loading vari- 
ance, with a red/green ratio observed on 
the cDNA microan-ay (Fig. 2a.c) for fi-actin 
ofl.CM. 
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UnlGcnc EST clustering sys- 
tem-'--. The scctind brj:e$i pnuip tif 
clones iCTiisisteJ o\ 1X3 scquenccil cDNA 
clone% j-eiwraied hy Mihiraciton cif cDNA 
lrt»m I he LhrumoMime*6 suppressed 
iuin-tumorti:eniL UACC-y03 l+ft) tell 
line vviih eDNA frnm its parental tumori- 
genie tell line UACtC-yo^ (rel. 
Appri>xim;iiely MM) additional genes 
(total 870 genes arrayed! were ohtained 
from KS r lihraries on the Iwsis «»f' I heir 
expre^sitin |utiern (I issue speeilk. and m) 
on), luieh array included the Uillowing 
hyhridi/aiion lonirols: plasmid vector, 
lamhtla. 0X 17-1 phage, loial human I )NA. 
human CaM\ DNA, and |>oly (AT. The 
syniheiic MaiulanU used lor normali/.)- 
lion ol' signals tn each wavelenglh were 
also arrayed. (!i»nirols were included in 
each t|nailram of the array to assess ihe reprnducihilitv til ihe 
hyhridi/ation signal. Two plaies old )NA cioiK*s dlerived I'rom 
the UAt;(;-*>ll3 suhiracied library! were alM» arrayed in dupli- 
cate, l-iileliiy ol'ihe Unigene array relative to dhK.S"r was tested 
l>y seijuencing ol' a random Minple of 1 1 cloiu-s ii.sed Itir 
microarray construe! ion. All seiiuentes were identical with the 



corres|>»nding dbl.S1 enineN. Addiiionalh. each inuitMi t.u. J 
cDNA Irnm ihe 1:AU:-*NI3 suhiracied lihrarv wa^ suiuen.cd. \ 
listing ol cDNAs compriMng ihis micniarrav uhuh vxcu 
diTivcd Irom the Unigene und hou.seUeping" panel v.m i>e 
ohuincd frnm hiip://wrtv.nih.guv/|)IU/LCr./AKKA^ cxpn - 
html. 

Northern blot analysis. Tola! RNA, 10 Mg |vr lane, wa^ elec- 
trophoresi'd in l.^%ag3n»e*lormatdchvde gels and transferred 
onto nylon membrane (H>4v»nd-N*, AmershamI by capillary 
blotting tnvrnighl. For DNA prolx's insert fragments from the 
Soares INIU cDNA library'" were obtained by vfcior PCR for 
p2l, MAUCK.S, a-l-aniichymoirypsin and (J-aciin. Hrobes lor 
fibroblast collagenase and MCAH/MllH- 1 were tsolated from a 
UACC-90.M+(^! enriched cDNA library*' with all pr.»bes 
labelled by random priming. Fillers were washed to a siiin- 
genc>' tif O.lx SSC at A2 *C lor 2U min. 

Websites. htip://cmgm.sianlord.edu/pbrown for pmiocols and 
pans list pertaining to micrturray fabrication. 
hlip://www.nchgr.nih.gt»v/l)|R/|.(:c;/AUI{AY/expn.html for a 
listing of cDNAs comprising this microarray which were 
derived from the Unigene ami 'housekeeping panel. 
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